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TECHNICAL  ACCOMPLISHMENTS 


A.  INTRODUCTION 

In  general,  the  research  program  has  continued  in  the  directions  outlined  in  the 
original  proposal  except  that,  because  of  the  departure  of  Dr.  Y.  G.  Yu  in  February,  six 
months  earlier  than  expected,  activity  in  the  development  of  practical  theoretical  methods 
for  the  calculation  of  the  conductivity  and  nonlinear  optical  properties  of  polymers  was 
severely  curtailed  during  the  last  six  months  of  the  project  however  two  manuscripts  have 
been  written  on  this  topic  (Appendices  II  and  III).  However,  considerable  progress  has 
been  made  in  obtaining  results  from  accurate  calculations  using  ab  initio  methods  (papers 
1-3,  Appendix  I). 

B.  TECHNICAL  PROGRESS 

1.  Ab  Initio  Calculations  of  Nonlinear  Optical  Properties  for  Small  Molecules. 

Although  there  are  a  large  number  of  semi  empirical  calculations  on  nonlinear 
optical  parameters  available  in  the  literature  there  have  been  until  very  recently  very  few 
accurate  ab  initio  calculations  of  any  kind  on  even  diatomic  molecules.  Our  philosophy 
is  that,  since  the  measurement  of  nonlinear  optical  properties  per  se  is  so  difficult,  the 
only  way  to  put  the  semi  empirical  methods  on  a  firm  foundation  is  to  provide  accurate 
reference  calculations  on  small  molecules  which  can  both  be  used  to  evaluate  the  results 
of  semi  empirical  methods  and  also  they  can  act  as  a  guide  in  choosing  parameters 
during  the  development  of  new  semi  empirical  methods  and  in  extending  the  present 
methods  to  the  heavier  elements.  This  task  was  regarded  as  one  of  the  main  aims  of  this 
research  contract.  Calculations  on  the  folloviring  sets  of  molecules  have  been  completed 
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and  the  results  are  in  the  process  of  being  analyzed  fully  and  written  up  for  publication. 

1.  Ab  initio  calculations  on  the  water  molecule  Published.  Theochem.  234, 173-183 
(1991) (Appendix  1,  paper  3). 

2.  Ab  initio  calculations  on  hydrogen  cyanide  and  phosphaethyne  (H-CsP). 
Published.  J.  Phys.  Chem.  24,  5586-5589  (1990) (Appendix  1,  paper  2). 

3.  Ab  initio  calculations  on  ammonia,  formamide,  pyrrole  and  pyridine.  Manuscript 
in  preparation.  (A  preliminary  version  containing  the  results  is  attached,  Appendix 
II). 

4.  Ab  initio  calculations  on  molecules  of  the  type  HgX  =  Y  p(  =  C,  Si,  Ge;  Y  =  O,  S). 
Manuscript  in  preparation. 

5.  Ab  initio  calculations  on  phosphine  and  arsine  (NHg,  PHg).  In  preparation. 

6.  Ab  initio  calculations  on  molecules  of  the  type  HgX  p<  =  S,  Se,  Te).  In  preparation. 

7.  Ab  initio  calculations  on  molecules  of  the  type  XOg  (X  =  S,  Se,  Te).  In 
preparation. 

8.  Ab  initio  calculations  of  model  molecules  containing  silicon  and  germanium,  e.g., 
of  the  type  H-XsX-H  and  XgH^.  In  preparation. 

9.  Ab  initio  calculations  of  small  molecules  containing  transition  metals.  Preliminary 
results  only.  Work  still  in  progress. 

These  calculations  were  designed  primarily  to  give  answers  to  the  following 
questions: 

1.  How  important  are  d  and  p  polarization  functions  on  the  hydrogens? 

2.  How  important  is  the  role  of  diffuse  p  functions  on  atoms  other  than  hydrogen? 
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3.  How  small  a  basis  set  can  we  get  away  with  without  degrading  the  values  of 
nonlinear  parameters? 

4.  Can  we  approximate  the  core  electrons  by  an  effective  potential  without  losing 
accuracy  and  thereby  save  CPU  time  and  reduce  memory  requirements? 

5.  Are  the  basis  sets  found  to  be  successful  for  small  molecules  transferable  to  larger 
molecules? 

6.  How  important  is  electron  correlation? 

7.  Are  there  any  general  rules  which  can  be  deduced  form  our  experience  which  will 
help  in  future  calculations? 

8.  How  do  the  accurate  results  obtained  compare  with  results  of  the  semi  empirical 

methods  now  in  use?  ’  - 

Calculations  on  the  hydrides  of  groups  three,  four,  five  and  six  fully  confirm 
our  prediction  that  there  is  a  dramatic  increase  in  the  value  of  the  second 
hyperpolarizability  as  we  ascend  a  particular  row  of  the  periodic  table.  The 
calculations  on  compounds  containing  the  heavier  elements  confirm  our  earlier 
results  on  extraordinary  sensitivity  of  the  /3  and  y  values  for  the  heavier  elements 
to  the  presence  of  very  diffuse  basis  functions  on  both  the  central  atom  gnd  the 
hydrogens.  For  example^n  calculations  on  compounds  of  galliunr^it  was  found  that 
one  extra  set  of  diffuse  p  functions  over  those  we  thought  would  be  adequate 
changed  the  value  of  y  by  a  factor  of  four. 

The  most  difficult  problem  we  have  had  to  solve  is  the  calculation  of  accurate 

values  of  all  the  first  hyperpolarizabilities  0?)  in  one  calculation.  Our  results  for  the 
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hydrides  (e.g.,  the  paper  3  in  Appendix  I  on  the  water  molecule)  show  that  the  p  values 
are  very  sensitive  to  the  number  and  type  of  polarization  functions  placed  on  the 
hydrogens.  In  particular^t  was  found  necessary  for  the  heavier  elements  to  have  multiple 
sets  of  both  p  and  d  functions(very  diffuse  function^) present.  Moreover^ these  functions 
seem  to  interact  in  a  complicated  way  so  that  erratic  values  of  certain  p  values  are 
obtained  by  leaving  off  one  or  more  functions.  Also  it  was  found  that  this  sensitivity  was 
different  for  different  elements.  For  example^the  calculation  on  the  p  parameters  for 
selenium  compounds  require  a  different  set  of  polarization  functions  on  the  hydrogen  than 
for  arsenic.  Reliable  values  can  be  obtained  by  saturating  the  hydrogens  with  five  sets 
each  of  d  and  p  functions  with  exponents  of  1.0,  0.5,  0.25,  0.125,  0.0625  but  for  a 
molecule  like  Ai^Hg  even  at  the  Hartree-Fock  level  this  becomes  prohibitive  in  terms  of 
computer  time  and  disk  space.  Moreover^  with  multiple  sets  of  polarization  functions 
converge  problems  occur.  These  problems  make  it  impossible  to  perform  this  type  of 
calculation  on  a  main  frame  computer  which  is  used  by  others^since  the  calculations  are 
I/O  bond^and  have  to  reside  in  the  computer  for  more  than  seven  days,  which  means 
they  are  often  lost  due  to  system  crashes  or  the  taking  over  of  the  computer  by  the 
sytsem  programmers  for  maintenance  and  updating.  We  have  overcome  this  difficulty  to 
some  extent  by  purchasing  a  IBM  RISC/6000  work  station  and  dedicating  it  for  long 
periods  of  time  to  the  calculation  of  nonlinear  optical  properties.  However  our  limited 
budget  prevented  us  from  purchasing  an  external  disk  drive  larger  than  1 .2  gigabytes  so 
we  are  somewhat  restricted  in  the  size  of  the  basis  set  we  can  employ  (<200  basis 
functions).  However,  despite  these  problems^ we  were  able  to  obtain  fairly  accurate 
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values  for  the  molecules  listed  above.  This  is  a  tremendous  advance  over  the  values  in 
the  literature  which  are  restricted  to  formaldehyde  (CHgO),  water  (HgO),  ammonia  (NHg) 
and  few  diatomics  and  most  of  these  have  appeared  in  the  last  year  or  so. 

If  we  can  be  satisfied  with  less  accurate  /3  values  (say  with  errors  of  30%-50%  'i'o  r 
may  be  one  matrix  element  with  the  others  within  20%)  then  we  have  found  that  reliable 
values  of  the  dipole  polarizabilities  (a)  and  hyperpolarizabilties  (y)  can  be  obtained  fairly 
easily  for  molecules  containing  more  than  a  dozen  atoms.  This  is  exemplified  by  our 
calculations  on  pyridine  and  pyrrole  in  the  manuscript  in  Appendix  II.  We  are  now 
pursuing  this  approach  vigorously  for  molecules  containing  the  heavier  elements.  The 
problem  here  is  that  we  need  realistic  geometries.  However,  we  believe  we  have  solved 
this  problem  as  described  later  in  this  report. 

Anotehr  aspect  of  our  work^  which  we  have  been  the  first  to  investigate 
systematicall^is  the  role  of  electron  correlation.  Sadlej  (e.g.,  Theoretica  ChimicaActa  79, 
123-160  (1991))  demonstrated  the  importance  of  including  electron  correlation  for  dipole 
polarizabilities.  His  group  has  been  using  many  body  perturbation  theory  (MBPT)  to 
fourth  order.  We  have  applied  Moller-Plesser  perturbation  theory  to  second  order  (MP2) 
and  have  obtained  very  similar  results  except  that  we  have  also  got  data,  for  the  first  time, 
for  the  effect  of  including  electron  correlation  for  the  first  hyperpolarizabilties  ()9)  and  we 
have  started  work  on  the  second  hyperpolarizabilties  (y).  The  results  for  /3  show  that  the 
effect  of  including  correlation  is  much  greater  for  these  parameters  than  it  is  for  the  a 
values. 

Indeed,  one  of  our  main  conclusions  is  that  it  is  a  very  difficult  and  time  consuming 
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task  to  get  accurate  values  for  even  the  smallest  molecule.  However  it  is  possible  to  get 
approximate  values  and  this  should  be  pursued  in  preference  to  rushing  ahead  to  get 
values,  which  are,  probably  meaningless,  fftm  a  semi  empirical  method  which  has  been 
inadequately  tested  out  on  smaller  molecules. 

Since  we  consider  the  production  of  accurate  values  of  nonlinear  optical 
parameters  of  reference  molecules  to  be  a  necessary  requirement  before  the 
development  of  semi  empirical  methods,  our  attention  in  the  latter  half  of  1991  was 
focused  on  whether  reliable  ab  initio  values  of  nonlinear  optical  parameters  could  be 
attained  using  smaller  basis  sets  than  we  found  necessary  for  very  accurate  calculations. 
One  method  we  have  explored  is  to  describe  the  behavior  of  the  inner  electrons  by  an 
effective  core  potential  (ECP  method).  Progress  in  this  area  is  described  in  the  next 
section. 

2.  Effective  Core  Potential  Calculations 

There  are  many  reports  in  the  literature  on  methods  for  the  development  of 
effective  core  potentials  (ECPs)  as  a  substitute  for  the  explicit  treatment  of  the  inner 
electrons  (i.e.,  those  with  a  different  (lower)  principal  quantum  number  than  that  of  the 
valence  electrons)  which  are  thought  to  be  responsible  for  the  major  part  of  the  binding 
in  a  molecule.  However,  there  have  been  relatively  few  attempts  to  actually  use  the  ECP 
method  for  the  calculation  of  molecular  properties,  other  than  the  energy.  Since  our 
interest  is  in  the  calculation  of  electric  field  effects  and  these  are  sensitive  primarily  to  the 
outermost  electronic  density  it  seemed  reasonable  to  try  out  the  use  of  ECPs  for  the 
purpose  of  carrying  out  calculations  of  nonlinear  optical  parameters. 
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Although  the  study  is  still  in  its  infancy  enough  results  have  been  obtained  to 
indicate  that  this  is  a  viable  way  to  proceed,  especially  for  molecules  containing  heavier 
element^  since  the  saving  in  computer  time  and  disk  space  becomes  relatively  greater 
the  greater  the  ratio  of  the  number  of  core  electrons  to  valence  electrons  becomes.  At 
the  same  time  as  we  have  investigated  the  use  of  an  ECP  we  have  also  attempted  the 
dangerous,  but  nevertheless  necessary,  procedure  of  reducing  the  number  of  diffuse 
functions  used  in  the  calculations  of  larger  molecules^  by  using  our  experience  with  the 
accurate  calculations  to  choose  only  those  functions  found  to  be  most  sensitive  to  the 
applied  field.  The  results  given  in  Appendix  II  for  formamide  (HCONHg)  are  typical  of  the 
results  obtained.  Formamide  was  chosen  because  a)  it  is  a  very  polar  molecule  with  a 
large  amount  of  charge  transfer  between  a  donor  (NHg)  group  and  an  acceptor  (C  =  0) 
group  which  should,  if  qualitative  ideas  on  the  relationship  between  structure  and  large 
hyperpolarizabilities  have  any  validity,  result  in  enhanced  y  and  p  values  compared  to  the 
values  of  the  parent  molecules  (ammonia,  formaldehyde);  b)  there  are  previous  semi 
empirical  calculations;  and  c)  there  is  considerable  interest  in  having  a,  p,  and  y  values 
for  this  molecules  as  a  prototype  building  block  for  polypeptides.  These  results  show  that 
using  the  ECP  method  results  in  a  considerable  saving  in  computer  time  with  little  change 
in  the  numerical  accuracy  of  the  results.  Similar  results  have  been  obtained  for  larger 
nitrogen  compounds  (e.g.,  pyrrole  and  pyridine)(see  Appendix  II). 

3.  Beyond  the  Hartree-Fock  Zero  Frequency  Static  Field  Methodology 

Before  one  can  expect  a  good  one  to  one  relationship  between  experiment  and 
theory  for  nonlinear  optical  parameters,  one  must  have  calculations  which,  in  addition  to 
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an  adequate  basis  set,  must  also  include  electron  correlation  and  also  which  include  the 
dynamic  frequency  dependence.  Our  computer  codes  have  been  modified  to  include 
both  these  effects  although,  at  the  moment,  calculations  on  the  effect  of  electron 
correlation  on  the  second  hyperpolarizability  have  to  be  obtained  by  performing  a  number 
of  separate  calculations  containing  electric  fields  of  varying  strength  at  a  number  of 
different  positions  and  then  using  the  appropriate  expansion  formulae,  a  very  laborious 
process.  Typical  of  the  results  obtained  are  those  we  obtained  for  -fo  r  Uy  i.J  -f* 
(Appendix  II,  Table  2)  for  including  electron  correlation  at  the  MP2  level.  While  expecting 
an  increase  in  computer  time  we  were  somewhat  unprepared  for  the  large  increase 
experienced  even  for  Ck-n  atom.  It  was  also  found  that  it  was  necessary  to  increase  the 
convergence  criterion  on  the  total  energy  by  one  or  two  orders  of  magnitude  to  get 
numerically  accurate  values  of  the  first  polarizability.  (This  finding  leads  us  to  question 
on  the  grounds  of  numerical  accuracy  alone  many  literature  p  values  obtained  using  semi 
empirical  methods,  since  these  usually  employ  higher  tolerances  for  energy  convergence 
than  we  have  found  to  be  necessary.  The  values  obtained  in  these  calculations,  and 
others,  show  that  including  the  electron  correlation  changes  the  hyperpolarizabilities  by 
a  factor  between  1.5  and  2.0. 

4.  Geometry  Determination 

In  calculations  on  the  more  common  molecules^such  as  formamide  or  pyridine^ 
experimental  geometries  are  available  and  can  be  used,  but  for  the  more  exotic  molecules 
containing  tellurium  or  arsenic  experimental  geometry  studies  are  relatively  rare.  In 
principle  one  can  perform  theoretical  geometry  optimizations  using  either  ab  initio  or  semi 
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empirical  techniques.  However,  semi  empirical  techniques  for  molecules  containing 
heavier  atoms  are  notoriously  unreliable.  Also  there  have  been  few  ab  initio  studies  on 
molecules  containing  more  than  five  or  six  atoms  and  on  molecules  containing  atoms  of 
the  third  and  higher  periods  (e.g.,  our  study  on  the  geometry  of  SO2  (Appendix  I,  paper 
1).  Moreover,  there  have  been  very  few  attempts  to  use  ECP  gradient  methods  to  obtain 
geometries.  For  these  reasons  we  have  been  looking  at  methods  for  the  calculation  of 
geometries.  As  soon  as  we  find  more  that  6-10  atoms  in  a  molecule  we  ran  into 
computational  problems.  This  was  traced  to  a  lack  of  stability  in  the  initial  Hassian  matrix 
used  to  calculate  the  force  constants  needed  to  start  the  geometry  optimizations 
especially  for  cyclic  molecules.  A  new  method  was  developed  and  implemented  by  Dr. 
Yu.  This  required  an  initial  semi  empirical  (AM1  method)  geometry'  optimization,  followed 
by  a  force  constant  calculation.  The  force  constants  were  then  used  as  the  starting  point 
for  the  ab  initio  calculation.  Even  though  the  calculated  semi  empirical  geometry  could 
be  poor  the  force  constants  calculated  using  it  were  found  to  facilitate  the  convergence 
of  the  ab  initio  geometry  optimization  considerably  and  cut  down  the  computer  time 
required  by  one  or  two  orders  of  magnitude.  This  has  enabled  us  to  perform  geometry 
optimizations  on  relatively  large  molecules  and  as  benzothiophen  CgHgS.  Examples  of 
applications  of  this  procedure  are  being  written  up  separately  from  the  nonlinear  optical 
parameter  studies,  but  an  adequate  geometry  optimization  method  is  a  necessary 
prerequisite  before  the  hyperpolarizability  calculations  are  even  begun. 

Another  problem  we  encountered  in  the  geometry  calculations  was  the  problem 
of  what  basis  set  to  use.  Even  though  the  basis  set  requirement  of  geometry 
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optimizations  is  nowhere  near  the  stringent  requirements  of  the  nonlinear  optical 
parameter  calculations,  some  polarization  functions  are  needed.  There  have  however 
been  few  studies  made  on  how  to  choose  these  functions.  We  have  therefore  extended 
our  study  on  the  SOg  molecule  to  study  molecules  of  the  general  type  MO^,  where  M  is 
S,  Se,  Te,  Cr  and  Mo  and  n  is  1,  2  or  3.  As  a  check  on  the  calculations  we  have  also 
performed  vibration  frequency  calculations  since  matrix  isolation  experiments  of  these 
quantities  are  often  more  available  for  "exotic"  molecules  than  are  geometry  studies. 
These  results  have  been  generally  good  (e.g.,  TeOg)  however  the  results  on  the  geometry 
of  SeOg  have  been  uniformly  poor.  In  every  calculation  we  have  tried,  including  a  very 
extensive  configuration  interaction  calculation,  the  bond  distance  has  always  been  much 
shorter  than  experiment  even  though  the  vibration  frequencies  have  been  reasonable. 
This  has  lead  us  to  suspect  the  experimental  values  rather  than  the  theoretical  methods 
so  we  are  proceeding  to  other  molecules,  leaving  this  molecule  alone  until  more 
experience  has  been  obtained  on  the  type  of  basis  set  we  have  to  employ  to  get 
reasonable  geometries.  Calculations  on  other  small  molecules  containing  transition  metal 
atoms  are  either  completed  and  awaiting  analysis  or  are  in  progress. 

5.  Ab  initio  calculation  of  Polymers 

The  other  main  thrust  of  the  research  program  was  to  develop  a  computer  code 
for  the  calculation  of  the  electronic  structure  of  organic  polymers.  The  method  chosen 
was  one  developed  over  the  last  two  decades  by  the  research  group  of  Professor  J.  J. 
Ladik  at  the  Friederich  Alexander  University  at  Erlangen  -  Nurnberg,  Germany.  This 
program  has  been  used  previously  to  investigate  vibrational  and  transport  properties  of 
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organic  polymers  and  is  the  only  program  which  incorporates  electron  correlation  into  its 
formalism.  We  have  installed  this  program  at  the  University  of  Tennessee  and  Professor 
J.  Ladik  visited  us  for  a  period  during  November  and  December  (1989)  to  instruct  us  in 
matters  concerning  the  use  of  the  program.  We  have  used  this  program  to  calculate 
band  structures  of  polymers  containing  elements  from  the  third  and  higher  groups  of  the 
periodic  table.  As  might  be  expected,  for  a  problem  of  this  complexity,  the  calculations 
are  not  easy  to  perform  and  it  is  very  easy  to  get  false  results.  The  program  for  example 
is  very  sensitive  to  the  geometry  and  is  very  difficult  to  get  converged  results  with  the 
larger  basis  sets  we  would  like  to  use!  In  addition  to  the  problem  of  basis  set  size  we 
also  have  the  problem  of  how  many  interacting  neighbors  are  to  be  included.  Typical  of 
the  results  we  have  been  able  to  complete  is  a  study  of  the  electronic  band  structure  of 
polyethylene,  polysilene  and  polygermane  using  two  basis  sets,  i.e.,  a  minimal  basis  and 
a  double  zeta  set.  A  manuscript  describing  this  work  is  in  preparation  (see  Appendices 
III  and  IV).  The  main  feature  of  the  work  is  that  it  demonstrates  the  dramatic  drop  in  band 
gap  in  going  from  the  carbon  polymer  (24.0, 15.6  eV)  to  the  silicon  polymer  (14.1, 10.2 
eV)  and  the  germanium  polymer  (9.3  eV,*).  (The  minimum  basis  result  is  given  first,  the 
double  zeta  value  second.)  (*  means  that  for  the  Ge  polymer  indicates  we  were  not  able 
to  get  convergence  after  many  tries!)  These  values  cannot  be  compared  with  experiment 
quantitatively  because  we  do  not  yet  have  correlation  effects  in  our  method.  However  the 
trend  should  be  qualitatively  correct.  Analysis  of  the  orbitals  involved  in  the  band  gap 
show  that  although  Sc  and  Ge  behave  similarly  they  are  considerably  different  from 
carbon.  The  results  support  the  idea  that  in  these  molecules  there  is  considerable 
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delocalization  of  the  sigma  electrons. 

6.  The  Calculation  of  Nonlinear  Properties  of  Polymers 

No  one  has  yet  succeeded  in  developing  a  working  method  for  the  calculation  of 
the  nonlinear  optical  properties  of  polymers.  Ladik  recently  suggested  a  formalism  for 
doing  these  calculations  and  since  we  already  had  the  necessary  polymer  code  the 
extension  of  the  method  in  this  direction  was  a  reasonable  project  to  attempt.  This  was 
a  very  formidable  task  but  was  completed  by  the  end  of  the  summer  (1990)  by  my 
colleague,  Dr.  J.  G.  Yu.  Since  then  calculations  have  been  completed  and  compared  with 
calculations  on  small  similar  molecules  calculated  by  the  Hartree  Fock  methods  described 
in  Appendix  III.  The  basic  difficulty  in  the  methods  proposed  for  calculating  electric  field 
effects  of  polymers  is  that  because  of  the  term  |e|  E^Z  (E2  being  the  electric  field  in  the 
2  direction)  the  chain  will  become  non-periodic.  In  order  to  get  around  this  difficulty 
Ladik  suggested  that  the  polymer  chain  of  2N^  +  1  unit  cells  be  divided  into  2Ng  +  1 
segments  (Ng  <  N^).  Each  segment  is  then  allowed  to  contain  a  certain  number  of  unit 
cells  each  of  which  are  given  an  average  value  for  the  dipole  moment  component  Zf^. 
The  success  of  the  method  then  depends  upon  the  insensitivity  of  the  parameter  to  the 
way  the  partitioning  was  made.  Our  initial  calculations  on  polyethylene  and  its  silicon  and 
germanium  analogues  using  a  minimal  basis  set  were  quite  successful  (Appendix  III)  but 
we  are  running  into  problems  when  we  increase  the  size  of  the  basis  set  and  at  the 
present  time  we  have  not  come  up  with  a  satisfactory  scaling  method  which  is  completely 
independent  of  basis  set. 

The  computer  code  has  been  extended  to  calculate  the  second  hyperpolarizabilities 
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(y)  and  has  also  been  written  and  tried  out  using  minimal  basis  sets  purely  as  a  means 
of  checking  the  numerical  accuracy.  As  expected  for  such  a  small  basis  set  the  results 
are  physically  meaningless  so  they  are  not  quoted  here. 
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INTRODUCTION 

There  is  currently  a  great  interest  in  the  design  of  organic  materials  with  certain 
characteristic  nonlinear  optical  properties.  In  order  to  understand  the  relationship 
between  the  electronic  structure  of  molecules  and  their  nonlinear  properties  a  number  of 
different  semi  empirical  methods  have  been  proposed.  In  order  to  place  these  methods 
on  a  sounder  footing  and  also  to  extend  their  parameterization  to  elements  for  which  very 
little  experimental  information  is  available  we  are  systematically  calculating  accurate  values 
of  the  dipole  polarizabilities  and  the  first  and  second  hyperpolarizabilities  for  small 
molecules  using  the  ab  initio  finite  field  method  is  implemented  in  the  HONDO  and  the 
Gaussian  computer  codes.  The  calculations  reported  here  are  for  a  number  of  nitrogen 
compounds,  namely  ammonia,  formamide,  pyridine  and  pyrrole.  For  ammonia  the  results 
of  very  large  basis  sets  are  reported  and  compared  to  some  recent  similar  calculations 
using  slightly  smaller  basis  sets.  It  is  then  shown  that  very  similar  results  can  also  be 
obtained  using  smaller  basis  sets.  The  smaller  basis  sets  are  then  used  for  the  larger 
molecules  and  the  results  are  compared  with  the  results  of  semi  empirical  calculations. 
RESULTS  AND  DISCUSSION 
A.  Formamide 

Formamide  was  chosen  as  an  example  of  a  medium  sized  molecule  because  there 
is  widespread  interest  in  the  polarizabilities  and  hyperpolarizabiiities  of  this  molecule  as 
a  starting  point  for  the  development  of  parameters  for  the  studies  on  intermolecular 
interactions  of  proteins  and  DNA.  The  molecule  is  also  of  interest  as  an  example  of  the 
enhancement  of  nonlinear  optical  parameters  due  to  charge  transfer,  i.e.,  between  NH2 
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as  an  electron  donor  and  the  carbonyl  group  as  an  electron  acceptor. 

The  results  of  using  different  basis  sets  are  reported  in  Table  1.  All  calculations 
were  performed  at  the  HF  6-31 G  optimized  geometry.  The  reference  calculation  I  (column 
2)  was  performed  at  the  HF  level  with  a  3-21 G  basis  set  supplemented  on  the  C,  0  and 
N  by  one  set  of  diffuse  s,  two  sets  of  p  and  3  sets  of  d  polarization  functions  using  the 
exponents  suggested  by  Dykstra  et  al.  A  similar  calculation  was  also  performed  in  which 
a  second  diffuse  p  set  (exp  =  0.005)  was  included.  As  in  the  case  of  ammonia  this 
omission  had  negligible  effect  and  the  results  given  in  Table  1  are  for  a  calculation  which 
does  not  include  this  diffuse  p  function.  The  second  p  set  was  also  omitted  in  the  other 
calculations  reported  in  Table  1.. 

Our  ultimate  aim  is  to  perform  calculations  on  larger  molecules  than  formamide  so 
the  effect  of  replacing  the  core  electrons  by  an  effective  core  was  investigated.  These 
results  are  given  in  column  3  (Calc  II).  Except  for  the  very  sensitive  and  py^ 
elements  there  was  surprisingly  little  change  but  the  computer  time  required  was 
considerably  reduced  (by  30%).  The  ECP  method  was  also  used  for  Calc  III  and  IV.  In 
Calc  III,  Calc  II  was  supplemented  by  two  sets  of  p  functions  on  the  hydrogen.  Although 
this  addition  doubled  the  computation  time  it  had  very  little  effect  on  the  dipole 
polarizability  tensor.  There  were  however  considerable  changes  in  several  of  the  p 
values,  expecially  the  p^  term  which  was  changed  by  34%. 

In  Calc  IV  a  d  set  (exp  =  0.1)  was  added  to  the  hydrogens.  There  was  a 
negligible  change  in  a  and  only  relatively  small  changes  in  the  first  and  second 
hyperpolarizabilities,  indicating  that  acceptable  results  can  be  obtained  without  going  to 
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the  expense  of  adding  the  d  set.  In  Calc  V  a  very  small  set  of  polarization  functions,  i.e., 
only  one  set  each  of  d  and  p  functions  with  exponents  of  0.05  as  suggested  by  DuPuis 
et  al.  for  larger  molecules  was  used.  There  was  a  loss  in  5  of  only  11%,  however,  there 
were  large  changes  in  the  p  values  showing  that  the  results  for  the  first  hyperpolarizability 
using  small  basis  sets  could  be  unreliable,  y  also  changed  by  more  than  about  20%, 
indicating  that  first  and  second  hyperpolarizabilities  for  polar  molecules  using  a  limited 
number  of  polarization  functions  may  only  have  qualitative  significance. 

Finally,  in  order  to  assess  the  effort  of  including  electron  correlations  the 
polarizability  and  first  hyperpolarizability  (HF  and  MP2)  calculations  were  performed  using 
the  G90  computer  code.  The  basis  set  was  the  standard  3-21 G  valence  set 
supplemented  by  the  polarization  functions  for  carbon  and  nitrogen  and  hydrogen 
recommended  by  Dykstra  (a  total  of  120  contracted  functions).  The  results  are 
summarized  in  Table  2.  The  HF  calculation  is  very  similar  as  expected  to  the  ECP 
calculation  III  of  Table  1.  It  can  be  seen  that  the  effect  of  including  electron  correlation 
at  the  MP2  level  is  considerable  on  all  the  properties.  The  reduction  of  the  dipole  moment 
by  10.5%  brings  it  into  excellent  agreement  (3.71  Debyes).  The  average  dipole 
polarizability  is  also  increased  by  12.7%,  again  bringing  it  into  excellent  agreement  with 
an  experimental  value  quoted  by  Waite  et  al.  (J.  Chem.  Phys.  83,  4047  (1985).  There  are 
no  experimental  results  for  the  hyperpolarizabilities  with  which  to  compare  results  for  the 
)9  elements.  This  is  unfortunate  as  the  changes  in  these  native  elements  including 
electron  correlation  are  quite  dramatic.  However  one  must  conclude  that  it  is  very 
important  to  either  include  electron  correlation  explicitly  or  to  make  a  correction  for  it 
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before  comparing  with  experiment. 

B.  Pyridine 

In  order  to  assess  the  effect  of  destroying  the  symmetry  of  the  benzene  molecule 
on  the  nonlinear  optical  parameters,  especially  the  js  values  which  are  course  zero  by 
symmetry  in  benzene,  a  number  of  calculations  were  performed  on  pyridine.  These  are 
compared  with  calculations  on  benzene  in  Table  3.  Columns  2  and  3  are  the  results  of 
a  benzene  (1/9)  and  a  pyridine  (IP)  calculation  performed  with  the  same  basis  set  namely 
the  standard  3-21 G  valence  basis  on  the  carbon,  nitrogen  and  hydrogen  supplemented 
by  one  set  of  p  and  one  set  of  d  functions  on  the  carbon  and  nitrogens  (a  similar  set  has 
been  used  for  recent  literature  calculations  on  benzene  (J.  Phys.  Chem.  95,  5873  (1991) 
and  other  enjinged  hydrocarbons  (J.  Phys.  Chem.  91,  4728  (1989):  J.  Chem.  Phys.  91, 
4728  (1989)).  The  results  predict  a  small  decrease  In  the  average  dipole  polarizability 
largely  due  to  a  reduction  in  the  out  of  plane  matrix  element  and  a  drop  in  the  average 
second  hyperpolarizability  (y)  of  about  20%.  Column  4  contains  the  results  (IIP)  of  adding 
a  set  of  p  functions  (exp  =  0.1)  to  the  hydrogens.  There  is  a  small  increase  in  the  dipole 
polarizability  and  a  further  drop  of  6%  In  the  second  hyperpolarizability.  The  first 
hyperpoiarizabilities,  as  expected,  are  more  sensitive,  but  even  there  the  change  in  the 
most  sensitive  element  0^  is  only  13%.  In  column  5  are  the  results  of  a  much  larger 
calculation  (HIP)  in  which  three  sets  of  d  polarization  functions  were  used  on  the  carbon 
and  nitrogen  using  exponents  taken  from  the  triple  zeta  ELP  set  of  Dykstra  (i.e.,  0.9, 0.13 
and  0.02)  and  two  sets  of  p  polarization  functions  on  the  hydrogens.  This  made  a  total 
of  247  uncontracted  Gaussian  functions  contracted  to  207.  This  is  close  to  the  maximum 
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we  can  handle  and  consumed  over  three  hours  of  CPU  time  on  our  IBM  3090/400 
supercomputer  to  get  the  polarizabilities  and  first  hyperpolarizabilities.  (Since  the  G86 
code  was  used  for  this  calculation  we  were  unable  to  get  the  second 
hyperpolarizabilities.)  As  expected,  there  are  no  dramatic  changes  in  the  results  for  the 
dipole  polarizabilities  but  all  three  /3  elements  were  increased  in  magnitude  with  the  largest 
change  being  in  0^,  which  was  40%  larger  than  the  result  of  Calc  1P.  The  results  of  all 
the  calculations  of  the  dipole  polarizabilities  are  also  in  good  agreement  with  the 
experimental  values  listed  in  column  7.  This  is  in  sharp  contrast  to  the  results  of  a  recent 
semi  empirical  calculation  (column  8).  This  calculation  which  is  typical  of  most  semi 
empirical  calculations  in  the  literature  which  almost  universally  tend  to  grossly 
underestimate  the  out  of  plane  contributions  to  a  because  they  do  not  contain  the  proper 
diffuse  p  and  d  functions.  Since  our  calculations  have  shown  that  the  diffuse  function 
requirements  for  the  calculation  of  0  and  y  values  are  even  more  demanding  with  respect 
to  the  number  and  type  of  functions  added  the  semi  empirical  methods  currently  in  use 
will  yield  very  unreliable  results  for  these  parameters  whereas  even  the  small  basis  set  we 
have  used  seems  to  give  semi  quantitatively  correct  results. 

C.  Pyrrole 

Pyrrole  is  a  molecule  of  great  interest  with  respect  to  its  electronic  properties 
because  it  is  contained  in  many  polymers  with  interesting  semiconducting  and  nonlinear 
optical  properties.  There  have  been  no  previous  ab  initio  calculations  on  the  nonlinear 
optical  properties.  Our  results  using  three  different  basis  sets  are  summarized  in  Table 
4.  The  results  for  a  and  y  are  very  similar  to  tiiose  obtained  for  benzene.  Certainly,  there 
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is  no  enhancement  shown  by  putting  in  the  NH  group  instead  of  a  carbon-carbon  bond. 
The  results  show  that  there  is  little  to  be  gained  by  including  the  p  functions  with  exp  = 
0.9  on  the  hydrogens  but  that  there  is  considerable  change  in  the  p  values  by  using  1  set 
of  p  functions  (exp  =  0.1)  on  the  hydrogens.  This  is  unfortunate  since  their  addition 
increases  the  size  of  the  basis  set  by  1  function.  Thus  if  values  of  p  are  needed  (i.e.,  with 
an  error  less  than  50%)  a  molecule  the  size  of  pyrrole  is  probably  the  limit  that  is  practical 
at  the  present  time,  especially  if  as  is  important  electron  correlation  is  included. 
CONCLUSIONS 

The  major  conclusions  of  this  work  are  that  reasonable  ab  initio  results,  especially 
if  an  effective  core  potential  (ECP)  is  employed,  can  be  obtained  for  the  dipole 
polarizabilities  and  the  second  hyperpolarizabilities  of  quite  large  molecules  using  a  split 
valence  (3-21 G)  basis  set  supplemented  by  two  diffuse  functions  (d  and  p).  On  the  other 
hand  much  larger  sets  are  needed  to  get  values  of  the  second  hyperpolarizabilities 
accurate  to  within  fifty  per  cent  even  at  the  Hartree  Fock  level.  The  results  on  formamide 
also  show  that  electron  correlation  is  extremely  important  for  the  first  hyperpolarizabilities, 
much  greater  than  for  the  dipole  polarizabilities.  The  results  also  show  that  current  semi 
empirical  methods  are  not  even  qualitatively  correct  especially  for  the  out  of  plane 
contributions.  (Note:  In  the  final  manuscript  a  table  of  comparison  with  results  obtained 
using  the  new  MOPAC  semi  empirical  method  wll  be  included.) 
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Table  1. 


Nonlinear  Parameters  (a,  /3,  y)  for  Formamide. 
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HI® 
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V® 
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4.207 

4.194 

4.190 

4.169 

4.55 
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25.510 

25.729 

25.900 

26.010 
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18.037 

18.410 

18.509 

18.413 
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28.55 

28.798 

28.980 

29.090 

29.080 

d 

24.033 

24.311 

24.463 

24.505 

21.588 

Aa 

^xxx 

39.605 

41.404 

41.633 

45.310 

60.792 

^XX2 

17.507 

18.800 

19.633 

21.035 

28.310 

^XYY 

5.536 

7.404 

9.908 

9.088 

12.720 
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7.949 

8.071 

7.646 
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8.762 

^Yft 

17.777 

15.821 

14.827 

15.103 

20.783 

^777 

-19.131 

-18.913 

-20.143 

-19.750 

-39.525 

Y 

3175.86 

3257.23 

3343.98 

3455.50 

5514.68 

®HF  standard  32-1 G  set  with  3D  sets  of  polarization  functions  on  C  and  N  (Dykstra 
exponents),  no  pol  functions  on  H  (102  basis  functions). 

^ECP  calculation  same  valence  set  as  Calc  I. 

®ECP  calculation  as  in  Calc  II  but  with  2  p  pol.  functions  on  H  (exp  =  0.9,  0.1). 

^ECP  calculation  as  in  Calc  III  but  with  1  d  set  added  to  H  (exp  =  0.1). 

®ECP  calculation  using  1  p  and  1  d  set  (exp  =  0.05)  on  the  C  and  N. 


Table  2. 


MP2  Calculations  on  the  Nonlinear  Parameters  of  Formamide. 


HF® 

MP2® 

4.190 

3.750 

25.870 

28.928 

2.88 

3.456 

18.190 

20.033 

28.852 

33.195 

24.304 

27.385 

44.110 

66.740 

20.687 

23.170 

8.761 

0.122 

7.426 

12.179 

15.180 

28.546 

-20.477 

-16.476 

^Standard  3-21 G  basis  set  supplemented  by  2  p  sets  on  H  (exp  =  0.9,  0.1 
respectively) (120  basis  functions). 


Table  3.  HFSCF  Nonlinear  Optical  Parameters  (a,  p,  y)  for  Benzene  and  Pyridine. 


IB® 

IP® 

IIP^’ 

HIP® 

IVP'^ 

VP® 

VIP‘ 

0.00 

2.80 

2.79 

2.40 

“xx 

73.83 

71.28 

72.35 

74.82 

75.46 

76.9 

67.34 

“yy 

73.83 

66.097 

67.07 

69.55 

70.07 

70.2 

59.85 

“zz 

38.398 

34.67 

36.54 

40.72 

40.83 

37.1 

8.43 

a 

62.02 

57.35 

58.65 

61.69 

62.13 

61.4 

45.2 

^xxy 

0.00 

20.34 

18.54 

21.42 

21.04 

^zzy 

0.00 

16.80 

16.90 

21.32 

18.16 

^yyy 

0.00 

-10.347 

8.98 

-14.53 

-10.40 

Y 

14,595.7 

11,415.0 

10,731 

®Basis  set  I;  HF3-21G  with  1  set  of  p  and  1  set  of  d  functions  (exp  =  0.050)  on  the 
carbon  and  nitrogen  atoms. 

‘’Basis  set  II:  Basis  set  I  and  1  set  of  p  functions  (exp  =  0.1)  on  each  hydrogen,  total  of 
133  functions. 

®Basis  set  III:  HF3-21G  with  1  diffuse  s  and  1  diffuse  p  and  a  3d  set  on  carbon  and 
nitrogen  as  in  the  ELP  set  of  Dykstra  et  al. 

‘‘Basis  set  IV:  Basis  set  III  supplemented  by  two  sets  of  p  functions  on  hydrogen  (exp 
=  0.9,  0.1)  as  for  the  ELP  set. 

®Experimental  results  quoted  in  Mol.  Physics  28,  577  (1979). 

Values  are  recent  semi  empirical  values  calculated  using  the  MNDO  method  {Theochem. 
22S.  193  (1991). 


Table  4.  Nonlinear  Optical  Parameters  tor  Pyrrole. 


11*^ 

lll*^ 

IM  (Deb) 

1.917 

1.858 

1.867 

57.33 

60.91 

61.97 

“yy 

33.606 

38.44 

39.00 

®77 

56.035 

59.26 

60.43 

a 

48.991 

52.87 

53.80 

^xxz 

-28.88 

-37.76 

-40.244 

^yyz 

17.310 

33.58 

34.704 

^777 

20.99 

35.101 

36.33 

Y)oo(x 

21,322.58 

18,260.97 

Yioyy 

9,864.56 

10,951.3 

4,257.78 

3,331.2 

YyYYY 

13,030.68 

17,968.4 

Yyyzz 

5,859.08 

5,722.7 

Y7777 

10,007.97 

7.345.6 

Y  (AV) 

16,864.82 

16,717.1 

®HF  3-21 G  with  1  set  of  p  and  1  set  of  d  polarization  functions  on  C  and  N  (exp  = 
0.05) (1CX)  basis  functions). 

^HF  as  in  Calc  I  but  with  3  sets  of  d  functions  (exp  =  0.9,  0.03,  0.02)  on  C  and  N  and  1 
set  of  p  (exp  =  0.1)  on  H  using  the  ELP  set  of  diffuse  functions  of  Dykstra  (180  basis 
functions). 

‘^HF  (G.86)  as  in  ^  but  with  2  sets  of  functions  (exp  =  0.9,  0.1)  and  1  set  of  s  functions 
(exp  =  0.01)  on  H  (200  basis  functions). 
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ABSTRACT 


The  approach  suggested  by  Ladik  [6]  for  calculating  polarizabilities  of  polymers  is 
implemented  to  calculate  the  polarizabilities  of  polyethylene,  polysilane  and  polygermane 
using  a  minimal  basis  set.  The  results  using  the  same  method  of  segmentation  for  all 
three  polymers  are  in  satisfactory  agreement  with  accurate  ab  initio  results  for  the 
corresponding  oligomers. 
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1. 


Introduction 


Because  polymers  are  now  being  recognized  as  having  potential  as  a  class  of 
materials  for  electronics  and  optical  applications,  there  has  been  considerable  research 
in  quite  different  scientific  areas  in  an  attempt  to  design  polymers  with  the  described 
properties.  The  design  of  polymers  with  suitable  non-linear  optical  properties  is  one  area 
which  has  recently  been  given  a  good  deal  of  attention  by  both  theoreticians  and 
experimentalists. 

In  the  dipole  approximation,  the  microscopic  response  of  a  material  to  a  uniform 
static  electric  fie’d  E  can  be  described  as 

Pi  (E)  =  i^i  *  ^  •  •  (1) 

where  M|  is  a  component  of  the  induced  molecular  dipole  moment,  a  is  the  linear 
polarizability  tensor  and  p  and  y  are  quadratic  and  cubic  hyperpolarizability  tensors, 
respectively. 

Within  the  one-electron  restricted  Hartree-Fock  (RHF)  approximation,  two  basically 
different  approaches  can  be  used  to  calculate  the  static  (hyper)polarizabilities  of 
molecules  and  oligomers:  the  finite-field  (FF)  and  the  sum-over-state  (SOS)  methods^  . 
Basically,  the  SOS  and  FF  methods  correspond  to  different  approximative  levels  of  the 
perturbed  Hartree-Fock  theory.  There  have  also  been  attempts  to  calculate 
(hyper)polarizability  tensors  on  polyene  oligomers  as  models  for  long  chain  polymers^’^ 
using  ab  initio  methods.  These  calculations  however  are  very  expensive  and  are  limited 
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to  the  number  of  momomer  units  that  can  be  utilized,  e.g.,  15  monomeric  units  in  the 
case  of  polyenes^. 

Quite  recently,  Barbier  et  al"*  tried  to  use  SOS  scheme  to  calculate  the 
polarizabilities  of  polymers  directly.  They  adopted  a  perturbative  method  previously 
developed  by  Gerkin  and  Mednis®.  In  their  calculations  on  a  model  molecular  hydrogen 
chain  and  also  on  polyacetylene  they  used  STO-3G  basis  set  and  found  that  the 
polarizability  along  the  axis  (a^^)  of  the  polymers  were  similar  to  those  from  calculations 
on  the  corresponding  oligomers. 

The  FF  approach  formally  equivalent  to  the  coupled  perturbed  Hartree-Fock 
(CPHF)  theory  is  basically  superior  to  SOS  method  because  it  allows  for  the  electrons  to 
relax  self-consistently  in  the  presence  of  the  perturbing  electric  field.  There  is  however 
a  fundamental  difficulty,  defined  previously  by  Ladik  in  the  IF  approach.  This  is  taht  the 
addition  of  the  external  electric  and  one  is  not  troubled  by  orbiting  decisions  on  the  type 
and  number  of  excited  states  to  include  in  the  direction  of  the  polymer  chain  destroys  the 
periodic  symmetry  of  the  chain,  which,  unfortunately,  is  a  requirement  of  the  rigorous 
application  of  ab  initio  Hartree-Fock  Crystal  Orbital  (CO)  theory  to  polymers.  Ladik® 
suggested  an  approximate  method  to  solve  the  problem,  in  which  the  neighbor  unit  cells 
of  the  quasi-dimensional  polymers  were  combined  as  segments  and  within  each  segment 
the  values  of  the  electric  field  are  taken  as  a  constant,  but  no  applications  were  reported. 

In  the  present  paper,  Ladik’s  method  is  implemented  for  the  first  time  and  tested 
by  calculating  the  polarizability  of  some  polymers. 

2.  Theory 
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According  to  the  ab  initio  SCF  LCAO  CO  (Hartree-Fock)  theory  of  periodic  chains 
using  Born-von  Karman  boundary  conditions^’®  for  a  linear  polymer  with  z  being  the 
polymer  axis,  it  is  necessary  to  solve  the  following  generalized  matrix  eigenvalue  equation: 

F(k)  Ci(k)  -  ii(k)S{k)Ci(k)  (2) 


with 


Slk)  -  ^  e^KgaS(Q):  «  <x^|xJ> 


(3) 


and 


F(k)  -  f;  e^xgaFig)  ;  [F(g)]„  =  <x;|i?*'  +  £  |xS> 

r-K  pi 


(4) 


In  the  presence  of  a  homogenerous  electric  field  E,  the  additional  term  |  e  |  rE  will 
occur  in  (eqn.  (4)).  The  new  Fock  matrix  is  then 

F^k.m  -  y;  e^F(g.E)  AFkg.E)}^,  -  <Xr\fi’’  *  |e|r^  + 

«•-»  (5) 

^  -  [F(<7)]„  ♦  |e|<x;|(rr|x?> 

{XXXXXXXXXXXXXXXX},  Ladik®  has  proved  that  the  supplemental  term 
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(6) 


E  e^***<Xrl  |e|  zE  \xl> 


remains  finite  if  a  Gaussian  basis  is  used  despite  the  fact  that  v  =  |e|  rE  is  unbounded 
and  the  presence  of  E  does  not  destroy  the  hermiticity  of  the  Fock  matrix: 

F(y..E)*  *  F(x.E) 

A  more  difficult  problem  introduced  by  the  existence  of  lejrE  is  the  aperiodicity 
caused  by  E^.  For  an  aperiodic  chain,  the  whole  formulas  given  by  eqns  (2)  -  (4)  are  no 
longer  valid.  Ladik  suggested  an  approximate  treatment  for  the  aperiodicity  to  surmount 
the  difficulty.  The  suggested  procedure  is  to  assume 

|e|^.E  -  |e|^,<Xrkix;>  £  i^sinCJcga)  <Xrklxl> 


for  the  case  that  E^  is  not  very  large,  instead  of  calculating  the  matrix  elements  <Xr‘’ 
I  EjjZ  I  •  if  ca'^  done  dividing  the  chain  of  2N  +  1  unit  cells  into  2N  +  1  segments 
(N  <  N)  where  each  segment  has  a  certain  number  of  unit  cells.  In  each  segment  m,  z 
has  an  average  value  Z^.  If  a  constant  number  of  cells  is  used  for  the  segment  length 
t,  then  2N  +  1  =  «(2N  +  1)  and  eqn.  (8)  becomes 

Therefore,  in  the  presence  of  an  electric  field 
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®"'®*<Xr|2|xl>  “  l®|£^,<Xrklxl>  +  J^2i8in(Jcqra)<x;irlxl>  *  le|£,[fs,5(o) 

q—K  fl-1  '  ' 

H 

^  2isin  (kqa)  ] 

g-1 


iFik,  £)„•  £  e"®'<X/|^'  +  T  *  \e\{E^  *  E^)  lx?>  + 

q~-H  px 

H  .  B 

|e|r,t2iS'„(0)  +5'2i  sin(Jcga)  4  ] 

n  "  ■»! 


(10) 


The  generalized  eigenvalue  equation  is: 

Fj  (k,  E)  Cj,  (k,  E)  -  Ci  Ik.  E)  Slk)  (k.  E)  (11) 

The  one-electron  wavefunctions  in  the  presence  of  an  electric  field  can  be  written  as 

N  M 

^iik.zE)  -  —l—Y^  e^*«*]rCi,,{Jr,£:)X?s<l>r(^^ 

V2W+i  pi 

where  m  is  the  number  of  basis  functions  per  unit  cell.  The  dipole  moment  per  unit  cell 
induced  by  an  electric  field  can  be  calculated  as 

jr 

«  n* 

^ind  .  |P|4>x(£)>  =  2^/  dkg  <(J>i(/c,£)  |er|<|)j(ic,£)>  (13) 

~ 

Sometimes,  another  definition  of  dipole  moment  can  be  given 
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(14) 


|1^  m  aE  +  1/2PE^  +  l/2lyE^E  +. .  . 


t-i  ^  * 


Comparing  above  both  equations,  the  tensor  elements  for  the  (hyper)polarizabilities  can 
be  calculated  as 


,lad 


(15) 


“  dEjdEj“0 


(16) 


ind 


“  dEjSsi^Ej  ■ 


(17) 


3.  Computer  Code 

From  the  viewpoint  of  programing,  in  order  to  calculate  polarizability  and 
hyperpolarizability  tensors  of  polymers,  an  additional  one-electron  intergral 

5^  e^*««<Xr||e|r£:|x|>  (18) 


is  added  to  original  one-electrc^^rjn^tri^._|cin.^1^  is^Ji|  sum  of  dipole  integrals 
weighted  by  the  appropriate  coefficients.  The  subroutines  calculating  dipole  integrals 
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were  supplied.  With  those  integrals  the  first  two  additional  integrals  introduced  by  the 
presence  of  electrical  field  are  obtained.  The  third  term,  which  is  the  most  important  for 
the  study  of  non  linear  optical  properties  of  polymer,  is  more  difficult  to  calculate.  It  is 
necessary  to  choose  the  parameters  2^,  and  N,  followed  by  the  summation  over 
overlap  integrals  given  in  eqn.  (10). 

With  newly  joined  Fock  matrix  F(k,  E)  and  solving  the  eigenvali  'iiquation  in  the 
presence  of  electrical  field  (eqn.  (11)), 

Fi{k.E)Ci{k,E)  ^  Ci(k,E)Slk)Ci(k.E)  (11) 

Sometimes,  with  the  help  of  dipole  integrals,  the  induced  dipole  moment  can  be 
easily  calculated. 

Numerical  differentiation  is  then  applied  to  get  polarization  and  hyperpolarization 
tensors.  The  usually  used  value  for  the  components  of  the  electric  field  is  0.001  a.u.  For 
example,  with  of  0.001  a.u.  can  be  calculated  as 

_  -i*i°** ( -g,)  ]  (20) 


The  calculation  of  and  Oyy  can  be  calculated  in  the  usual  way. 

4.  Results 

(a)  a^^andttyy 

Table  1  shows  comparative  results  of  and  calculated  with  the  polymer 
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program  with  HONDO  calculations  for  some  general  long-chain  oligomers. 

The  results  of  the  two  methods  for  all  three  systems  (polyethelene,  polysilane  and 
polygermane)  are  in  very  good  agreement.  However,  its  results  of  polymer  calculations 
are  all  smaller  than  the  HONDO  calculations.  The  and  Oyy  values  are  for  the 
H(C2H4)„H  unit  in  the  two  molecules,  n  equals  5  or  7.  Because  we  are  calculating  for  the 
infinite  linear  polymer  it  is  reasonable  to  assume  that  the  longer  the  model  molecule  used 
the  more  precise  the  result  obtained.  The  values  of  the  and  in  H(C2H4)7H 
calculation  are  less  than  those  of  H(C2H4)5H  and  the  results  on  the  polymer  have 
negative  errors  compared  to  those  of  H(C2H4)5H.  We  conclude  from  this  that  the 
approach  taken  does  produce  satisfactory  results  for  the  calculation  of  and  for 
polymers. 

(b)  The  choice  of  parameters  in  calculation  of 

A  much  more  difficult  problem  is  the  calculation  of  0^2,  the  polarizability  along  the 
axis  polymer  chain.  The  difficulty  is  that  the  presence  of  electric  field  this  direction 
destroys  the  periodicity  of  the  polymers  which  results  in  problems  in  solving  the  resulting 
Fock  equations.  Ladik’s  approximation  has  solved  the  problem.  In  implementing  the 
method  proposed  by  Ladik  it  is  necessary  to  choose  the  parameters  Z,  and  N  (m>  1), 
i.e.,  we  have  to  decide  how  to  divide  into  segments  and  how  to  choose  the  average 
electric  field  in  each  segment.  We  have  been  unable  to  arrive  at  a  theoretical  background 
for  choice,  so  that  we  have  had  to  resort  to  this  trial  and  error.  There  have  been  no 
previous  attempts  to  solve  this  problem  in  practice. 

Our  first  example  is  polyethylene.  The  contribution  for  of  both  diagonal  blocks 


PRELIMINARY  MANUSCRIPT 


of  electric  field  integral  <x°  |  |e|E2lXi‘’>  and  non-diagonal  blocks  <x°  |  lejE^IXs^^ 
have  been  analyzed  and  it  was  found  that  according  to  Ladik’s  approximation 

<Xrl|e|£,zlx>elel£,ZiS„(0) 


and  that  the  contribution  of  above  term  for  is  zero.  In  the  present  approahc  then  use 
of  2^(-1)'’Srs(0)  instead  of  <Xr“|e|E|z|Xs“>  when  N  is  even  then  <Xr“|2|Xs'’>  ^0. 
when  N  is  odd  then  <Xr“  |z|  Xs*’^  <  i-®-*  ^^®  of  the  integral  is  chosen  through  N. 

Table  3  presents  the  contributions  of  the  diagonal  and  non-diagonal  blocks  for 
for  polyethylene  using  a  STO-3G  basis  set.  In  these  calculations  we  used  1 .25  as  Z, , 
which  is  an  approximately  a  quarter  of  the  translation  vector  length.  These  results  show 
that  the  dominant  contribution  to  is  from  the  diagonal  blocks  while  the  non-diagonal 
blocks’  contribution  is  less  important.  In  these  above  calculations  three  neighbor  units 
were  taken  as  a  segment  and  the  choice  =  mZ.,  was  used.  Because  non-diagonal 
blocks  are  less  important  the  results  of  changing  the  method  of  dividing  the  segments 
have  no  remarkable  effect.  In  the  calculations  on  the  other  polymers,  we  adopt  the  same 
choice  as  for  polyethylene  for  the  needed  parameters. 

5.  Polarizabilities  of  polysilane  and  polygermane 

The  calculations  for  polyethylene  showed  that  a  quarter  of  the  translation  vector 
length  is  a  suitable  value  for  Zy  Using  this  choice  and  the  same  choice  for  Zy  Zm 
and  N,  of  polysilyane  and  polygermane  were  calculated.  The  suitable  values  of  Z  for 
polysilicane  and  polygermane  is  1.25  x  1/4  x  aQ.The  results  agree  well  with  the  accurate 
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calculations  on  the  oligor  .jrs. 


CONCLUSIONS 

The  results  are  the  first  attempt  to  calculate  the  polarizabilities  of  an  infinite  chain 
polymer.  The  results  for  polyethylene,  polysiline,  and  polygermane  seem  reasonable. 
They  show  the  enormous  increas  in  the  expected  as  carbon  is  replaced  by  silicon  and 
then  by  germanium. 

Attempts  are  in  progress  to  extend  this  approach  to  other  polymers  and  to  use  a 
more  realistic  basis  set. 
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Tabic  2. 


The  eflect  of  different  number  of  neighbors  in  calculations  of  a„  and  for  polyethylene 
oligomers. 


H-(QH,VH 

10.039 

9.941 

9.804 

9.690 

Table  3.  The  contributions  of  diagonal  and  non-diagonal  segments  to  a,*  in  the 
calculation  of  polyethylene. 


Number  of  Neighbors 

only  diagonal 

only  non- diagonal 

1 

-11.316 

1.121 

2 

-12.571 

1.115 

3 

-12.989 

1.122 

4 

-13.182 

1.124 

5 

-13.285 

1.125 

6 

-13.347 

1.125 

7 

-13.387 

1.125 

8 

-13.414 

1.125 

9 

-13.417 

1.132 

Table  4. 


Comparison  of  polarizabilities  of  polyethylene,  polysilane  and 
polygermane  with  polymers  and  oligomers  calculation. 


SiaH* 

GeaH*  | 

Oligomer 

Polymer 

Oligomer 

Pol)nmer 

Oligomer 

Polymer 

®xx 

10.039 

9.709 

19.920 

19.536 

25.023 

24.243 

“yy 

9.804 

9.368 

23.648 

23.145 

34.388 

33.220 

“ff 

12.829 

12.164 

46.298 

45.825 

109.747 

112.230 

Appendix  IV 


THEORETICAL  STUDY  ON  THE  ELECTRONIC  STRUCTURE  OF 
POLYETHYLENE,  POLYSILANE  AND  POLYGERMANE 

by 

Jianguo  Yu  and  John  E.  Bloor 


Department  of  Chemistry 
The  University  of  Tennessee 
Knoxville.  TN  37996-1600 


PREUMINARY  MANUSCRIPT 


Abstract 


The  electrical  structures  of  polyethylene,  polysilane  and  polygermane  have 
been  calculated  by  the  ab  initio  SCFCO  method  using  a  number  of  different  basis 
sets.  The  band  gaps,  band  widths,  ionization  potentials,  electron  affinities  and 
density  of  states  (DOS)  of  the  polymers  are  compared.  The  effects  of  different 
basis  sets  and  different  numbers  of  interacting  neighbors  on  the  energy  band 
figures  using  electrical  structures  of  the  polymers  are  discussed.  The 
conduction  mechanism  of  polysilane  and  polygermane  is  explored  using  energy  band 
figures . 


I .  Introduction 

There  is  a  very  great  interest  in  the  semiconducting  properties  of 
polysilanes  and  polygermanes  in  recent  years  because  of  the  potential 
technological  applications  of  these  materials.  The  understanding  of  these 
properties  requires  information  on  their  electronic  structure  -  especially  the 
band- edge  electronic  structure  and  the  conduction  band.  While  for  the  electrical 
structures  of  polyethylenes  and  polysilanes  there  have  been  previous  theoretical 
studies,  both  at  the  semi  empirical  (1-4]  and  ab  initio  [5-10]  level,  there  has 
been  for  polygermane  only  the  work  by  Takada  et.  al.  [1]  using  the  first 
principle  local  density  functional  method.  ^ 

In  this  paper,  the  ab  initio  RHFCO  method  developed  by  L^ik's  group  []  is 
used  in  the  comparative  study  of  electrical  structure  and  conduction  properties 
of  polyethylene,  polysilane  and  polygermane.  The  important  electronic  properties 
such  as  band  gap  and  band  widths,  ionization  potential,  electron  affinity  and 
density  of  states  (DOS)  of  those  polymers  are  compared.  The  effects  of  Including 
3d  orbitals  in  the  calculations  on  polysilane  is  discussed  and  the  Influence  of 
substituting  methyl  groups  into  polysilane  is  investigated.  The  conduction 
mechanism  of  polysilane  and  polygermane  is  also  explored. 

II.  Methods  and  Computational  Details 

The  energy  band  structures  of  the  one -dimensional  periodic  neutral  polymers 
were  calculated  using  the  ab  initio  Hartree-Fock  crystal  orbital  method  (12). 
To  obtain  information  on  how  the  electronic  and  energetic  properties  of  the 
macromolecules  depend  on  the  choice  of  the  basis  set  the  minimal  ST0-3G  [13], 
3-21G  [14],  6-31G  [15]  and  dementi's  double-zeta  [16]  basis  sets  were  used.  In 
some  calculations  on  Si  these  basis  sets  were  supplemented  by  d  polarization 


functions . 


For  all  the  three  polymers  the  planar  zigzag  structure  (figure  1)  was 


adopted. 


The  geometrical  parameters  of  the  polymers  used  in  the  calculations  are  given  in 
Table  1.  The  geometrical  parameters  of  polyethylene  are  taken  from  Otto  et,  al. 
(17].  The  optimized  values  for  the  geometry  of  polysilane  obtained  from 
molecular  mechanics  calculations  was  adopted  il8) .  For  polygermane  the  optimized 
values  from  ab  initio  3-21G  calculations  for  oligomers  of  H(Ge2H4)„H  (n  -  )  were 
used  to  construct  the  geometry  of  the  polymer  unit  cells  [19]. 

The  effect  of  methyl  substitution  for  silicon  polymer  was  studied  by 

\ 

comparing  the  electronic  structure  of  polysilane  and  poly(dimethylsilyne}  with 
different  basis  sets.  The  geometry  of  polydime thy Is ilyne  is  based  on  the 
molecular  mechanics  optimization  of  reference  18. 

All  calculations  were  performed  on  an  IBM  3090/200  computer  at  the 
University  of  Tennessee  Computing  Center  (UTCC)  using  the  ab  initio  polymer 
program  POLYMER  from  {  ) . 

III.  Results  and  Discussion 

1.  The  effect  of  the  different  number  of  interaction  neighbors. 

In  Table  2  the  theoretical  results  for  polyethylene  and  polys ilane 
in  terms  of  different  number  of  interacting  neighbors  with  both  ST0-3G  and 


■1 


dementi's  basis  sets  are  reported.  These  results  show  that  as  the  number 
of  neighbors  is  increased  for  both  polymers  (polyethylene  and  polysilane) 
and  with  the  difference  basis  sets  for  polyethylene,  the  valence  and 
conductive  band  gaps  are  shifted  toward  higher  energy.  However  the 
changes  in  the  width  of  these  bands  is  negligible.  The  changes  in  the 
ionization  potentials  and  electron  affinity  with  different  numbers  of 
interaction  neighbors  have  the  same  trends,  i.e.,  the  ionization 
potentials  become  smaller  and  the  electron  affinities  get  bigger.  The 
calculations  show  that  whereas  there  is  a  considerable  change  in  the  band 
gaps  in  increasing  the  number  of  nearest  neighbors  from  2  to  5  for  both 
polymers  and  for  both  basis  sets  for  polyethylene  it  is  the  further 
increase  from  five  to  eight  cause  little  further  change. 

2.  The  effects  of  basis  sets  on  electronic  band  structure. 

Table  3  presents  the  electrical  structure  of  polyethylene, 
polysilane  and  polygermane  with  different  basis  sets. 

For  polyethylene,  the  results  with  ST0-3G  and  dementi's  basis  sets 
show  big  differences.  The  valence  band  shift  to  lower  energy  about  1  eV 
and  its  width  becomes  narrower.  The  change  in  the  conduction  band  is  more 
marked.  It  is  lowered  by  about  11  eV  which  changes  the  band  gap  of  the 
polymer  from  23.97  eV  for  ST0-3G  to  14.60  eV  for  dementi's  basis  set. 
Most  of  this  change  is  in  the  electron  affinity  which  changes  from  14.07 
eV  to  3.29  eV. 

For  polysilane,  the  same  trend  as  polyethylene  is  exhibited.  In 
polysilane,  the  shift  of  valence  band  is  more  remarkable  (from  ST0-3G  to 
6-31G  it  shifts  toward  lower  energy  about  2.2  eV,  with  dementi's  basis 
set,  the  band  further  shifts  about  2.6  eV  toward  lower  energy).  The 
conductive  band  move  about  9  eV  toward  lower  energy  from  ST0-3G  to 
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dementi's  basis  sets  and  the  width  of  the  band  becomes  only  1.6  eV  from 
3.7  eV  of  ST0-3G.  The  more  noticeable  result  is  that  the  minimum  of 
conductive  band  with  dementi's  basis  set  becomes  negative,  which  makes 
the  electron  affinity  have  a  negative  value.  The  change  of  the  band  gap 
of  the  polymer  from  ST0-3G  to  dementi's  basis  set  isn't  as  big  as 
polyethylene  which  change  from  14.1  to  10.2  eV. 

Because  of  convergence  difficulties  with  the  large  basis  sets  only 
ST0-3G  and  3-21G  calculations  can  be  reported  for  polygermane.  The  same 
trend  of  the  electronic  structure  change  from  minimum  STO-3G  to  double - 
zeta  valence  3*21G  can  be  seen,  which  with  the  valence  3>21G,  the 
conductive  bands  shift  toward  lower  energy.  Because  the  basis  sets  only 
have  little  difference  the  change  of  electrical  structure  of  polygermane 
isn't  as  distinct  as  both  polyethylene  and  polysilane. 

3.  The  effect  of  3d  functions  in  the  calculation  of  electronic 
structure  of  polysilane  and  poly(dimethylsilylene) . 

For  the  system  Including  Si  atoms,  the  effect  of  3d  functions  still 

is  a  negliglbleless  problem.  For  polysilane  and  poly(dimethylsilylene) , 
\ 

the  energy  band-structures  were  calculated  with  STO-3G  and  ST0-3G*  (a  set 
of  3d  functions  of  Si  atom  with  exponent  -  0.39  was  supplied  to  ST0-3G 
basis  set)  and  the  results  are  presented  in  Table  4.  These  results  show 
that  the  addition  of  3d  function  on  the  Si  atoms  produces  little  effect  on 
the  band  structure  of  both  polysilane  and  poly(dimethylsilylene)  and  the 
effect  is  like  that  using  bigger  basis  sets  shown  in  last  section.  This 
shows  that  the  3d  functions  aren't  necessary  for  a  qualitative  comparative 
study  of  band  structure  of  polysilane  and  substituted  silicon  polymers. 

4.  The  comparison  of  polysilane  and  poly ( dimethyls ilylene ) . 

The  substitution  of  methyl  groups  into  polysilane  shifts  both  the 


valence  and  conductive  bands  about  1  eV  toward  higher  energy,  and  the 
widths  of  valence  bands  become  1  eV  narrower  while  the  conduction  band 
becomes  1.5  eV  wider.  The  band  gap  is  about  1  eV  narrower  in 
poly (dime thylsilylene)  than  in  polysilane. 

5.  A  comparison  of  conduction  properties  of  polyethylene,  polysilane 
and  polygermane. 

In  Table  3  are  listed  the  band  gaps  of  polyethylene,  polysilane  and 

polygermane.  The  experimental  value  of  band-gap  of  polyethylene  is  8.8  eV 

[20].  With  6-31G**  basis  set  and  both  RHF  and  HP  methods,  Sahai  [11] 

obtained  13.4  and  10.3  eV  for  the  band-gap  of  polyethylene  respectively. 

Our  result  with  dementi's  basis  set  is  very  similar  Sahai 's  RHF  value. 

It's  well -know  that  RHF  method  overestimate  the  band- gap.  In  our 

calculation  for  polysilane  with  dementi's  basis  set  the  band-gap  is  10.2 

eV.  Using  same  basis  set  and  RHF  CO  method  for  polyethylene,  calculated 

band-gap  is  about  6  eV  larger  than  experimental  that.  According  to  this 

proportion,  the  band-gap  should  be  about  4.2  eV,  which  is  approach  that 

gotten  by  other  ab  initio  and  serai  empirical  researchers  [1-10], 

\ 

For  polygermane  only  ST0-3G  and  3-21G  basis  sets  were  used.  The 
band  gaps  gotten  with  both  the  basis  sets  are  9.3  and  9.5  eV  respectively. 
In  terms  of  ST0-3G  basis  set,  the  band-gap  of  polygermane  is  4.3  eV 
smaller  than  that  of  polysilane.  It  can  be  expected  that  polygermane 
would  have  much  smaller  band  gap  than  polysilane.  In  the  local  density 
functional  calculation  [10]  by  Takada  et.  al.  the  difference  between  the 
band  gaps  of  polysilane  and  polygermane  is  only  0.6  eV. 

Despite  this  contradiction  we  can  still  draw  the  conclusion  that 
polysilane  and  polygermane  have  a  much  more  smaller  band  gap  than 
polyethylene,  although  the  atoms  C,  Si  and  Ge  of  those  polymers  are  in 


same  group  in  the  periodic  table.  The  difference  Is  the  essential 
property  which  makes  polysilane  and  polygermane  semiconductors  and 
polyethylene  a  typical  insulator. 

Figure  1  presents  the  energy-band  structures  of  polyethylene, 
polysilane  and  polygermane  with  ST0-3G  calculations.  This  figure  shows 
that  beside  the  difference  In  band-gaps  there  is  also  a  difference  in 
energy  band  profiles. 

For  each  polymer  the  first  six  highest  occupied  valence  bands  (HOVE) 
and  the  first  four  lowest  unoccupied  conduction  bands  (LUCE)  are  drawn. 
Solid  and  broken  lines  indicate  a  bands  and  x  bands,  respectively.  The 
HOVE'S  of  all  three  polymers  are  a  orbitals,  consisting  of  positive 
overlap  between  the  Pz  orbitals,  the  difference  between  polyethylene  and 
the  other  two  polymers,  polysilane  and  polygermane,  is  that  the  energy 
bands  of  polysilane  and  polygermane  are  highest  occupied  throughout,  while 
in  the  majority  of  K,  the  highest  occupied  valence  band  of  polyethylene  is 
a  K  band  which  consists  of  a  the  negative  overlaps  between  Pz  orbital  of 

both  carbon  atoms  and  between  Px  orbitals  of  each  carbon  atom  and  Is 

\ 

orbital  of  neighboring  H  atoms.  Nelson  and  Pietro  [21]  have  previously 
pointed  out  that  polysilanes  conduct  through  extensive  a  delocalization 
and  we  have  also  found  that  [19]  the  polygennanes  are  very  similar. 
Besides  having  bigger  band  gap,  the  energy  band  crossing  can  be  thought  as 
a  reason  polyethylene  is  an  insulator. 

From  an  analysis  of  the  coefficients  of  the  band  orbitals  it  can  be 
found  that  the  LUCE  of  polyethylene  are  different  in  composition  to  those 
of  both  polysilane  and  polygermane.  The  LUCB's  of  the  latter  two  consist 
of  negative  overlap  of  SP  hybridized  orbitals  of  the  Si  or  Ge  atoms,  while 
those  of  polyethylene  have  a  contribution  from  the  H  orbitals  besides  SP 


overlap . 

On  the  profile,  as  its  valence  band,  the  LUCB  of  polyethylene 
crosses  with  ir*  crystal  bands  while  those  of  polysilane  and  polygermane  do 
not.  All  of  above  those  can  be  explained  as  the  reason  that  polysilane 
and  polygermane  are  semi  conductors  while  polyethylene  is  a  insulator. 

6.  Density  of  states. 

The  density  of  electronic  states  (DOS)  for  polyethylene,  polysilane 
and  polygermane  are  shown  in  Figure  2.  The  calculation  of  the  DOS  only 
included  first  six  HOVB's  and  first  four  LUVB's  for  each  polymer.  The 
results  computed  with  the  help  of  the  ST0-3G  minimal  basis  set  have  been 
used  for  the  graphics.  It  can  be  seen  from  the  comparison  ol  the  DOS  of 
valence  band  that  those  of  polysilane  and  polygermane  have  similar  shape, 
although  the  position  of  peaks  of  polygermane  is  about  2  eV  lower  than 
those  of  polysilane.  The  shape  of  DOS  of  polyethylene  is  different  with 
those  of  polysilane  and  polygermane.  The  DOS  of  conduction  bands  of 
polyethylene  is  more  concentrated,  while  the  DOS  of  conduction  bands  of 

polysilane  has  a  maximum  about  12  Ev.  In  the  figure  of  the  DOS  of 

\ 

conduction  bands  of  polygermane  four  peaks  can  be  distinguished, 
corresponding  with  four  conduction  bands. 
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Table  1.  The  Geometrical  Parameters  (Distances  in  A  and  angles  in  degree) 
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<XXX 

<HXH 

(CjH*) 

1.54 

1.09 
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109.5 

(SiaH*) 

2.34 

1.48 
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107.4 
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2.47 

1.53 
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109.0 

Rc-h 

<HCH 

<SiSiC 

<SiCH 

(SiaCCHa),) 

2.35 

1.87 

1.10 

115.4 

108.5 

110.0 

Table  2.  The  effect  of  number  of  neighbors  (units  =  eV) 
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T^>ie  3.  The  band  characteristics  of  polyethylene,  polysilane  and  polygermaoe  using  different  basis  sets  (all  ener^  values  in  eV) 


Figure  1. 
Figure  2. 


Figure  Captions 

Band  structure  of  (a)  Polyethylene,  (b)  Polysilane  and  (c) 
Polygermane . 

Density  of  states  (DOS)  of  (a)  the  valence  band  and  (b)  the 
conduction  band  of  polygermane,  (c)  the  valence  band  and  (d)  the 
conduction  of  polysilane  and  (e)  the  valence  band  and  (f)  the 
conduction  band  of  polygermane.  The  horizontal  axis  is  energy  in 
eV. 


